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The transverse and longitudinal Doppler effects of the sunbeam spectra are shown to result in the earth 
parameters such as the earth-self rotation and revolution velocities, the earth orbit semi-major axis, the 
earth orbital angular momentum, the earth axial tilt, the earth orbit eccentricity, the local latitude and the 
mass of the sun. The sunbeam global positioning scheme is realized, including the earth orbital position. 
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Fig. 1 Sunbeam Global Positioning Device with a needle at 
the center of radial angle graph paper. 



Recently, I showed, using the sunbeam global 
positioning device shown in Fig. 1, the earth-self 
rotation and the earth orbiting around the sun 
and the novel macroscopic wave geometric 
effect of the sunbeam [1]. I obtained the earth 
axial tilt Q, the earth orbit eccentricity e and the 
earth orbital angle R s at winter solstice and the 
local latitude u, and made a calendar as well. 

However, the earth parameters such as the 
earth orbit semi-major axis a, the earth-self 
rotational velocity V E at the equator and the 
mass of the sun M s are assumed to be known. 
How about earth bound data of them? Yes, I 
show them. 



Here, I present a scheme to obtain the earth 
orbit revolution velocity, and the earth orbit 
angular momentum, a, V E R s , u, Q, e, and M s . 

The idea is a simple one, the Doppler shifts. 
The longitudinal (transverse) Doppler shift goes 
as linear (square) of the relative velocity 
between the observer and the signal source. Let 
me discuss the longitudinal Doppler shift first. 
In the morning (afternoon), the observer on the 
earth surface moves toward to (away from) the 
sun, resulting in the blue (red) shifts of the 
sunbeam spectra. At noon, the earth-self 
rotation does not yield any shift, since the earth 
rotation velocity is perpendicular to the 
sunbeam direction. The earth orbit radial 
directional velocity is finite, resulting in the 
longitudinal Doppler shifts. The transverse 
Doppler shifts are expected all the time. How 
large are the Doppler shifts? Numerical values 
are noted: one electron volt = 2.4 xl0 14 Hz. The 
earth-self rotation velocity at the equator is 
0.465 Km/sec = 1.55 xlO" 6 speed of light (C = 3 x 
10 5 Km/sec). The average earth orbital velocity is 
30 Km/sec = 10" 4 C and the earth orbit radial 
directional velocity of the order of 1.67 x 10" 6 C. 
The transverse Doppler shift of 1 electron volt 
light (1.24 micron) goes as 10" 8 x electron volt = 
2.4 xlO 6 Hz. The Doppler shifts of the sunbeam 
spectra are large enough to be observable by 
the present state of art of technology. 

To study the Doppler shift [2], I use the 
coordinate systems same as those used for the 



study of the sunbeam directional angle c 
velocity and recapitulate the results [1]. The bar 
on the earth surface is considered as the 
observer. I choose the earth center as the origin 
of coordinates as shown in Fig. 2. Let the earth- 
self rotating axis be the z-axis and the earth 
equatorial plane be the x-y plane. Let the z-y 
plane be the meridian, noon longitude, and b 
longitude measured from the meridian and u 
the latitude. Let the sunbeam from the sun 
with propagation unit vector come to the earth 

k = (0, - cosq, sinq) (Ik) 

making angle q with the y-axis. 

We consider here the spherical earth. The 
sunset or sunrise longitude g can be determined 
by the crossing longitude with the sunbeam 
front plane (sfp), when sfp hits the earth center 
making angle q with the z-axis. We consider the 
latitude u circle of radius cosu at z = sinu in the 
unit of the earth radius. The distance from the 
z-axis to sfp on the z-y plane is sinu tanq equal 
to cosu cosg, the projection of the radius cosu 
on the z-y plane as shown in Fig. 2. Thus, we 
obtain 

cosg = tanq tanu. (lg) 
This condition is for the radial unit vector 
r g = (cosu sing, cosu cosg, sinu) 

to be perpendicular to k, k • r g = 0. 

For a given u, the lengths of day and night are 
determined by the angle q by Eq. (lg). The g = 
and g = n indicate no day and no night, 
respectively. By Eq. (lg), the north pole region 
I u | > n/2 - | q | has no day (no night) for 
positive (negative) q, and the south pole region 
the other way around. At equinoxes, q = 0, g = 
n/2, we get equal day and night lengths. Strictly 
speaking, there is a small difference since q is 
changing from morning to sunset. At winter 
solstice, q is the earth axial tilt Q. The latitudes 
| u | = k/2 - I Q| are called Arctic and Antarctic 
circles, respectively. The latitudes u = |Q|are 



known as Tropic of Cancer and Capricorn, 
respectively. 
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Fig. 2 sfp: sunbeam front plane, krp: k-r plane with the 
sunbeam and a bar, rhp: r-h longitudinal plane, ahp: a-h 
plane, angle f between the sunbeam and the bar, angle c 
between krp and rhp, u: latitude, longitudes b at the bar 
and g at sunset or sunrise measured from the z-y plane, 
the meridian, S: length of shadow on ahp, B: length of the 
bar stands vertically on ahp. 

At the place where a bar stands, the radial 
unit vector r, the azimuthal unit vector a and 
the longitudinal unit vector h are given by 

r = (cosu sinb, cosu cosb, sinu), (lr) 

a = {- cosb, sinb, 0), and (la) 

h = (sinu sinb, sinu cosb, - cosu). (lh) 

We then get the shadow length S by a bar 
height B on a-h plane (ahp), with the angle f 
between the sunbeam and the bar, as 



S= Btanf , 



(2) 



cosf = - k • r = cosu cosq cosb - sinu sinq. (3) 

The shadow angle c between k-r plane (krp) and 
r-h longitudinal plane (rhp), is given by 

cose = a • (k x r) /sinf (4) 

= [sinu cosq cosb + cosu sinq]/sinf, or 
tanc = sinb/[sinu cosb + cosu tanq]. (5) 

For b = {0; g} = {Noon; Sunrise or Sunset}, we 
obtain the angles f and c as 

f = {u + q;7T/2}, (6) 

tanf = {tan(u + q); Infinite} , (7) 

c = {0; arc cos[sinq secu]}. (8) 

By Eqs. (2), (6) and (7), we can examine the 
earth spherical nature by measuring the 
shadow lengths at two locations on the same 
longitudinal line. We expect two different 
shadow lengths would show the spherical earth. 
We measure the earth radius r E with the arc- 
distance between two locations and the 
difference between latitudes. At sunrise or 
sunset, from Eq. (8), the angle c is 90 degrees 
only at equinoxes, q = 0. Sunrise and sunset 
directions are not always exactly east and west. 

Taking derivatives of angles f of Eq. (3) and c 
of Eq. (5) with b, we obtain for b = {0; g} = 
{Noon; Sunrise of Sunset}, 

df/dt = BTdf/db = {0 ; cosu cosq sing} BT, (9) 

dc/dt = BT3c/3b = {cosq csc(u+q); sinu} BT. (10) 

Here t is the sidereal time and BT = 3b/3t = 15 
deg/hr, the earth-self rotation angular velocity. 
There is no term of 3q/3b in Eq. (10) at noon. In 
other cases, it is negligibly small and neglected 
in Eq. (9) and Eq. (10). We estimate its size as 
following. The q changes 4x23.5 degrees in a 
year, dq/dt = 94/365.25 deg/day = 0.01072 
deg/hr and then dq/db = (5q/5t)/(5b/5t) < 
0.000715. Its omission is justified in the time 
frame of a daytime. 



The angle c velocity at sunset or sunrise (AcVS) 
by Eq. (10) is given as 

15 sin(local latitude) deg/sidereal hour. (US) 

This is same as the rotating rate of swing plane 
of Foucault pendulum [3], showing the earth- 
self rotation. The angle c velocity at noon (AcVN) 
by Eq. (10) is rewritten as 

15 [sinu + cosu cot{u+g)]deg/ hr. (UN) 

AcVN has an additional term resulted by the 
novel macroscopic geometric effect of the 
sunbeam making AcVN faster than AcVS [1]. 

At the local noon, the shadow has the 
shortest length and AcVN the maximum. The 
time difference between the local and standard 
time noons yields the local longitude = local 
standard longitude + BT x time difference hours. 
We have the sunbeam global positioning device. 

The earth orbital system around the sun 
is shown in Fig. 3. 1 realize the sunbeam 
direction in the earth system same as the 
radial direction in the earth orbital system. 

k = (0, - cosq, sinq) is r s = (cosR, sinR, 0) (12) 

in terms of the respective coordinate variables 
where the earth orbit angle R is measured from 
perihelion on the x-axis in the earth orbital 
system around S-sun. Moreover, the unit vector 
w of the earth-self rotation angular velocity BT 
in the earth system and its unit vector w s in the 
earth orbital system are given by, 

w = (0, 0, 1) and (13) 

w s = (sinQ cosR s , sinQ sinR s , cosQ), (14) 

where R s is the earth orbit angle R at winter 
solstice. Thus, we get easily [4] 

sinq = k«w = r s «w s = sinQ cos(R - R s ) or (15q) 



R = R s + arccos[sinq/sinQ]. (15R) 

The R s is obtained by the data of dq/dt at 
equinoxes by Eq. (19). 




\< n=a(i+B) * |-*-rp=a(i-e)-* 



Fig. 3 Earth orbit radial distance r(R) = r /(l + e cosR), R: 
the earth orbit angle, r = a(l - e 2 ), a and e are the semi- 
major axis and eccentricity of the ellipse, P: perihelion, VE: 
vernal (spring) equinox SS: summer solstice, A: aphelion, 
FE: fall equinox, WS: winter solstice, R s : angle at WS, Q: 
the earth axial tilt, W s = (sinQ cosR s , sinQ sinR s , cosQ): BT 
unit vector, r s = (cosR, sinR, 0): the r(R) unit vector. 

By Eq. (15q), I get the angle R velocity as 

dR/dt = -(dq/dt) cosq cscQ csc(R - R s ). (16q) 

This equals to that by the earth orbit [3] 

dR/dt = L/M E r 2 (R) 

= (L/M E r 2 )(l + ecosR) 2 , (16R) 

r(R) = r /(l + e cosR) , (16r) 

r = a(l-e 2 ) = L 2 /M E 2 GM S , (16a) 

L/M E r 2 = 27T:/Y(l-e 2 ) 3/2 , and (16L) 

Y = 27ra 3/ 7(GM S ) 1/2 . (16Y) 



the earth mass G the Newton gravitational 
constant, and r(R) the earth orbit radial distance 
from the sun, respectively. We consider the 
function Z(R) defined, for | R | < n/2, by 

Z(R) = [(dR/dt)(R)/(dR/dt)(R +n)] 1/2 

= r(R + :r)/r(R) (17r) 

= (l + ecosR)/(l-ecosR) (17R) 

= |(dq/dt)(R)/(dq/dt)(R + 7r)| 1/2 . (17q) 

Z(R) has a maximum value at R = 0, which can 
be found. By Eq. (17R) and Eq. (17q), we get 

e = [Z(0)-l]/[Z(0) + l]. (18) 

The eccentricity of the earth orbit is obtained 
with data of dq/dt at perihelion and aphelion. 
Moreover, by Eq. (15q) and Eq. (16q), dq/dt has 
local extreme values at equinoxes (q = 0), which 
can be found. For Z(R S + n/2) = Z E , we obtain 

R s = arcsin[(l-Z E )/e(Z E +l)]. (19) 

The angle R s at winter solstice is obtained with 
data of dq/dt at equinoxes. By integrating dR/dt 
of Eq. (16R), we obtain a calendar, D days from 
perihelion as 

D = (Y/2tt)[H - eE 1/2 sinR/(l + e cosR)], (20) 

H = 2arc tan[K tan(R/2)], 

E = l-e 2 , and 

K=[(l-e)/(l + e)] 1/2 . 

R and q are related by Eq. (15q) and Eq. (15R). 
The second term in Eq. (20) shows the earth 
elliptic orbit. 

The Doppler shift formula for frequencies Wj 
and W 2 in the systems 1 and 2 is given by [6] 



Here Y is the Kepler period of the elliptic orbit 
[5], L the earth orbital angular momentum, M E 



W 2 = gammaWi [1-beta] (21) 



beta = k«V 21 /C (21b) 

gamma = [1-{V 21 /C) 2 ] 1/2 , (21g) 

V 21 = (dr/dt)r s + (rdR/dt)v s + (V E cosu)a (21v) 

v s = e z x r s , r s x v s = e z = (0, 0, l) s , 

where V 2 i is the relative velocity between the 
systems 1 and 2. The unit vector e z is the z- 
direction in the earth orbital system, and its 
components (ei, e 2 , e 3 ) in the earth system are 
determined by the following conditions 

k»e z = =r s «e z , w»e z = e 3 = e z «w s = cosQ. (21z) 

The beta (gamma) results in the longitudinal 
(transverse) Doppler shift. 

Using Eq. (16R), Eq. (16r), Eq. (Ik) and Eq. (la), 
beta is obtained as 

beta = [VesinR - V E cosu sinb cosq]/C , (22) 

V=L/M E r , (23) 

V E = r E BT. (24) 

The first term in Eq. (22) is the earth orbit radial 
directional velocity and the second the earth- 
self rotation velocity projected on k, k»a V E cosu. 
At the local noon, b = 0, beta has only the first 
term. From perihelion to vernal equinox and to 
aphelion, we have the Doppler red shifts and 
the Doppler blue shifts in the rest region. 
First we consider the longitudinal Doppler shift 
and the sunbeam frequency F(P) at perihelion 
at noon (R=0, beta=0) as the reference 
frequency, disregarding the gamma factor. We 
write the frequency shift at R as 

W(R) = F(R) - F(P) = F(P)beta, (25W) 

W(VE) for at vernal equinox R=R VE =R s +7i:/2, etc.. 
Then, by Eq. (22) we get at noon, b = 0, 

Ve =C [W( Tt/2) - W(3n/2)]/2, (25V) 

tanR s = W(WS)/W(VE) = W(SS)/W(FE). (25R S ) 



Using Eq. (16a), Eq. (16L) and Eq. (23), we get 

a = V(l-e 2 ) 1/2 Y/2n, (27a) 

L = V 2 (l-e 2 ) 3/2 M E Y/2n, and (27L) 

M s = V 3 (1 - e 2 ) 3/2 Y/2nG. (27M S ) 

Using Eq. (lg), the second term in Eq. (22) may 
be rewritten at sunrise (SR) or sunset (SS) as 

-V E [cos(u+q)cos(u-q)] 1/2 /C. 

Neglecting small change in the first term in Eq. 
(22) during a daytime period, we get 

V E =C | W(SR) - W(SS) | (u=0,q=0)/2. (27V E ) 

The earth-self rotation angular velocity BT is 
obtained by BT = V E /r Ej with the earth radius r E . 
Using Eq. (25V) and Eq. (27V E ), we get 

[cos(u+q) cos(u-q)] 1/2 (any day) 

= C|W(SR)-W(SS)|/2V E (28) 

sinR =C[W(R at Noon) - W(P at Noon)]/Ve. (29) 

By Eq. (15q), Eq. (28) and Eq. (29), u and q are 
obtained. At winter solstice q = Q is obtained. 

We now discuss the transverse Doppler shift. 
For gamma of Eq. (21g), the relative velocity 
square is needed. For simplicity, we present 
the square of relative velocity at noon at 
perihelion as 

(V 21 ) 2 = V 2 (l+e) 2 + (V E cosu) 2 

-2VV E (l+e)cosu cosQ secq P , (30) 

sinq P = sinQ cosR s from Eq. (15q). 

The first term comes from the earth orbital 
revolution velocity and the second term from 
the earth rotational velocity. The last term 
comes from the scalar product of two velocities. 
If the sunbeam spectra in the rest frame are 



known by measuring He-lines or H-lines, we get 
the transverse Doppler shift at perihelion as 

W T (P) = I F(P) -F(rest) | = F(P)(V 21 /C) 2 /2 . (31) 

By Eq. (25V), Eq. (30) and Eq. (31), the earth 
orbit revolution velocity V and e are obtained. 
At aphelion, we get W T (A) by Eq. (30) with -e 
and q A =-qp- We get the frequency difference 

| F(A) - F(P) | =F(P)|(V 21 ) P 2 -(V 21 ) A 2 |/2C 2 

= F(P)2e(V/C) 2 [l-(V E /V)cosu cosQsecq P ]. (32) 

Here we used secq A = secq P by Eq. (15q). By Eq. 
(25V) and Eq. (32), V and e are obtained. 

I stress Eq. (15q) and Eq. (21z) play crucial 
roles for making connections between variables 
in the earth and orbital systems, and obtaining 
a calendar by q measured on the earth surface 
and the Doppler shifts of the sunbeam spectra. 

The Doppler shifts of the sunbeam spectra 
result in the parameters R s , V, V E and e without 
any parameter. One input variable Y is needed 
to get the earth orbit semi-major axis and the 
sun mass with G. For the earth orbit angular 
momentum, M E and Y are needed. Nowadays, 
M E is obtained by the satellite orbit semi-major 
axis and its period revolving around the earth. 

Moreover, the global positioning scheme is 
realized, including the earth orbit position on 
the local date, without any satellite. 

It is highly desirable to measure the Doppler 
shifts of the sunbeam spectra. Specially, the 
transverse Doppler shift of the sunbeam spectra 
would be physically one of challenging works. 

Photon has so many beautiful sides [7]. 
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